3.7 \(\int \frac{1}{(a+b \cot ^2(c+d x))^3} \, dx\)

Optimal. Leaf size=150 \[ \frac{\sqrt{b} \left (15 a^2-10 a b+3 b^2\right ) \tan ^{-1}\left (\frac{\sqrt{b} \cot (c+d x)}{\sqrt{a}}\right )}{8 a^{5/2} d (a-b)^3}+\frac{b (7 a-3 b) \cot (c+d x)}{8 a^2 d (a-b)^2 \left (a+b \cot ^2(c+d x)\right )}+\frac{b \cot (c+d x)}{4 a d (a-b) \left (a+b \cot ^2(c+d x)\right )^2}+\frac{x}{(a-b)^3} \]

[Out]

x/(a - b)^3 + (Sqrt[b]*(15*a^2 - 10*a*b + 3*b^2)*ArcTan[(Sqrt[b]*Cot[c + d*x])/Sqrt[a]])/(8*a^(5/2)*(a - b)^3*
d) + (b*Cot[c + d*x])/(4*a*(a - b)*d*(a + b*Cot[c + d*x]^2)^2) + ((7*a - 3*b)*b*Cot[c + d*x])/(8*a^2*(a - b)^2
*d*(a + b*Cot[c + d*x]^2))

________________________________________________________________________________________

Rubi [A]  time = 0.156941, antiderivative size = 150, normalized size of antiderivative = 1., number of steps used = 6, number of rules used = 6, integrand size = 14, \(\frac{\text{number of rules}}{\text{integrand size}}\) = 0.429, Rules used = {3661, 414, 527, 522, 203, 205} \[ \frac{\sqrt{b} \left (15 a^2-10 a b+3 b^2\right ) \tan ^{-1}\left (\frac{\sqrt{b} \cot (c+d x)}{\sqrt{a}}\right )}{8 a^{5/2} d (a-b)^3}+\frac{b (7 a-3 b) \cot (c+d x)}{8 a^2 d (a-b)^2 \left (a+b \cot ^2(c+d x)\right )}+\frac{b \cot (c+d x)}{4 a d (a-b) \left (a+b \cot ^2(c+d x)\right )^2}+\frac{x}{(a-b)^3} \]

Antiderivative was successfully verified.

[In]

Int[(a + b*Cot[c + d*x]^2)^(-3),x]

[Out]

x/(a - b)^3 + (Sqrt[b]*(15*a^2 - 10*a*b + 3*b^2)*ArcTan[(Sqrt[b]*Cot[c + d*x])/Sqrt[a]])/(8*a^(5/2)*(a - b)^3*
d) + (b*Cot[c + d*x])/(4*a*(a - b)*d*(a + b*Cot[c + d*x]^2)^2) + ((7*a - 3*b)*b*Cot[c + d*x])/(8*a^2*(a - b)^2
*d*(a + b*Cot[c + d*x]^2))

Rule 3661

Int[((a_) + (b_.)*((c_.)*tan[(e_.) + (f_.)*(x_)])^(n_))^(p_), x_Symbol] :> With[{ff = FreeFactors[Tan[e + f*x]
, x]}, Dist[(c*ff)/f, Subst[Int[(a + b*(ff*x)^n)^p/(c^2 + ff^2*x^2), x], x, (c*Tan[e + f*x])/ff], x]] /; FreeQ
[{a, b, c, e, f, n, p}, x] && (IntegersQ[n, p] || IGtQ[p, 0] || EqQ[n^2, 4] || EqQ[n^2, 16])

Rule 414

Int[((a_) + (b_.)*(x_)^(n_))^(p_)*((c_) + (d_.)*(x_)^(n_))^(q_), x_Symbol] :> -Simp[(b*x*(a + b*x^n)^(p + 1)*(
c + d*x^n)^(q + 1))/(a*n*(p + 1)*(b*c - a*d)), x] + Dist[1/(a*n*(p + 1)*(b*c - a*d)), Int[(a + b*x^n)^(p + 1)*
(c + d*x^n)^q*Simp[b*c + n*(p + 1)*(b*c - a*d) + d*b*(n*(p + q + 2) + 1)*x^n, x], x], x] /; FreeQ[{a, b, c, d,
 n, q}, x] && NeQ[b*c - a*d, 0] && LtQ[p, -1] &&  !( !IntegerQ[p] && IntegerQ[q] && LtQ[q, -1]) && IntBinomial
Q[a, b, c, d, n, p, q, x]

Rule 527

Int[((a_) + (b_.)*(x_)^(n_))^(p_)*((c_) + (d_.)*(x_)^(n_))^(q_.)*((e_) + (f_.)*(x_)^(n_)), x_Symbol] :> -Simp[
((b*e - a*f)*x*(a + b*x^n)^(p + 1)*(c + d*x^n)^(q + 1))/(a*n*(b*c - a*d)*(p + 1)), x] + Dist[1/(a*n*(b*c - a*d
)*(p + 1)), Int[(a + b*x^n)^(p + 1)*(c + d*x^n)^q*Simp[c*(b*e - a*f) + e*n*(b*c - a*d)*(p + 1) + d*(b*e - a*f)
*(n*(p + q + 2) + 1)*x^n, x], x], x] /; FreeQ[{a, b, c, d, e, f, n, q}, x] && LtQ[p, -1]

Rule 522

Int[((e_) + (f_.)*(x_)^(n_))/(((a_) + (b_.)*(x_)^(n_))*((c_) + (d_.)*(x_)^(n_))), x_Symbol] :> Dist[(b*e - a*f
)/(b*c - a*d), Int[1/(a + b*x^n), x], x] - Dist[(d*e - c*f)/(b*c - a*d), Int[1/(c + d*x^n), x], x] /; FreeQ[{a
, b, c, d, e, f, n}, x]

Rule 203

Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(1*ArcTan[(Rt[b, 2]*x)/Rt[a, 2]])/(Rt[a, 2]*Rt[b, 2]), x] /;
 FreeQ[{a, b}, x] && PosQ[a/b] && (GtQ[a, 0] || GtQ[b, 0])

Rule 205

Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(Rt[a/b, 2]*ArcTan[x/Rt[a/b, 2]])/a, x] /; FreeQ[{a, b}, x]
&& PosQ[a/b]

Rubi steps

\begin{align*} \int \frac{1}{\left (a+b \cot ^2(c+d x)\right )^3} \, dx &=-\frac{\operatorname{Subst}\left (\int \frac{1}{\left (1+x^2\right ) \left (a+b x^2\right )^3} \, dx,x,\cot (c+d x)\right )}{d}\\ &=\frac{b \cot (c+d x)}{4 a (a-b) d \left (a+b \cot ^2(c+d x)\right )^2}-\frac{\operatorname{Subst}\left (\int \frac{4 a-3 b-3 b x^2}{\left (1+x^2\right ) \left (a+b x^2\right )^2} \, dx,x,\cot (c+d x)\right )}{4 a (a-b) d}\\ &=\frac{b \cot (c+d x)}{4 a (a-b) d \left (a+b \cot ^2(c+d x)\right )^2}+\frac{(7 a-3 b) b \cot (c+d x)}{8 a^2 (a-b)^2 d \left (a+b \cot ^2(c+d x)\right )}-\frac{\operatorname{Subst}\left (\int \frac{8 a^2-7 a b+3 b^2-(7 a-3 b) b x^2}{\left (1+x^2\right ) \left (a+b x^2\right )} \, dx,x,\cot (c+d x)\right )}{8 a^2 (a-b)^2 d}\\ &=\frac{b \cot (c+d x)}{4 a (a-b) d \left (a+b \cot ^2(c+d x)\right )^2}+\frac{(7 a-3 b) b \cot (c+d x)}{8 a^2 (a-b)^2 d \left (a+b \cot ^2(c+d x)\right )}-\frac{\operatorname{Subst}\left (\int \frac{1}{1+x^2} \, dx,x,\cot (c+d x)\right )}{(a-b)^3 d}+\frac{\left (b \left (15 a^2-10 a b+3 b^2\right )\right ) \operatorname{Subst}\left (\int \frac{1}{a+b x^2} \, dx,x,\cot (c+d x)\right )}{8 a^2 (a-b)^3 d}\\ &=\frac{x}{(a-b)^3}+\frac{\sqrt{b} \left (15 a^2-10 a b+3 b^2\right ) \tan ^{-1}\left (\frac{\sqrt{b} \cot (c+d x)}{\sqrt{a}}\right )}{8 a^{5/2} (a-b)^3 d}+\frac{b \cot (c+d x)}{4 a (a-b) d \left (a+b \cot ^2(c+d x)\right )^2}+\frac{(7 a-3 b) b \cot (c+d x)}{8 a^2 (a-b)^2 d \left (a+b \cot ^2(c+d x)\right )}\\ \end{align*}

Mathematica [A]  time = 0.286408, size = 138, normalized size = 0.92 \[ \frac{\frac{\sqrt{b} \left (15 a^2-10 a b+3 b^2\right ) \tan ^{-1}\left (\frac{\sqrt{b} \cot (c+d x)}{\sqrt{a}}\right )}{a^{5/2}}+\frac{b (7 a-3 b) (a-b) \cot (c+d x)}{a^2 \left (a+b \cot ^2(c+d x)\right )}+\frac{2 b (a-b)^2 \cot (c+d x)}{a \left (a+b \cot ^2(c+d x)\right )^2}-8 \tan ^{-1}(\cot (c+d x))}{8 d (a-b)^3} \]

Antiderivative was successfully verified.

[In]

Integrate[(a + b*Cot[c + d*x]^2)^(-3),x]

[Out]

(-8*ArcTan[Cot[c + d*x]] + (Sqrt[b]*(15*a^2 - 10*a*b + 3*b^2)*ArcTan[(Sqrt[b]*Cot[c + d*x])/Sqrt[a]])/a^(5/2)
+ (2*(a - b)^2*b*Cot[c + d*x])/(a*(a + b*Cot[c + d*x]^2)^2) + ((7*a - 3*b)*(a - b)*b*Cot[c + d*x])/(a^2*(a + b
*Cot[c + d*x]^2)))/(8*(a - b)^3*d)

________________________________________________________________________________________

Maple [B]  time = 0.026, size = 363, normalized size = 2.4 \begin{align*} -{\frac{\pi }{2\,d \left ( a-b \right ) ^{3}}}+{\frac{{\rm arccot} \left (\cot \left ( dx+c \right ) \right )}{d \left ( a-b \right ) ^{3}}}+{\frac{7\,{b}^{2} \left ( \cot \left ( dx+c \right ) \right ) ^{3}}{8\,d \left ( a-b \right ) ^{3} \left ( a+b \left ( \cot \left ( dx+c \right ) \right ) ^{2} \right ) ^{2}}}-{\frac{5\, \left ( \cot \left ( dx+c \right ) \right ) ^{3}{b}^{3}}{4\,d \left ( a-b \right ) ^{3} \left ( a+b \left ( \cot \left ( dx+c \right ) \right ) ^{2} \right ) ^{2}a}}+{\frac{3\,{b}^{4} \left ( \cot \left ( dx+c \right ) \right ) ^{3}}{8\,d \left ( a-b \right ) ^{3} \left ( a+b \left ( \cot \left ( dx+c \right ) \right ) ^{2} \right ) ^{2}{a}^{2}}}+{\frac{9\,\cot \left ( dx+c \right ) ab}{8\,d \left ( a-b \right ) ^{3} \left ( a+b \left ( \cot \left ( dx+c \right ) \right ) ^{2} \right ) ^{2}}}-{\frac{7\,{b}^{2}\cot \left ( dx+c \right ) }{4\,d \left ( a-b \right ) ^{3} \left ( a+b \left ( \cot \left ( dx+c \right ) \right ) ^{2} \right ) ^{2}}}+{\frac{5\,{b}^{3}\cot \left ( dx+c \right ) }{8\,d \left ( a-b \right ) ^{3} \left ( a+b \left ( \cot \left ( dx+c \right ) \right ) ^{2} \right ) ^{2}a}}+{\frac{15\,b}{8\,d \left ( a-b \right ) ^{3}}\arctan \left ({b\cot \left ( dx+c \right ){\frac{1}{\sqrt{ab}}}} \right ){\frac{1}{\sqrt{ab}}}}-{\frac{5\,{b}^{2}}{4\,d \left ( a-b \right ) ^{3}a}\arctan \left ({b\cot \left ( dx+c \right ){\frac{1}{\sqrt{ab}}}} \right ){\frac{1}{\sqrt{ab}}}}+{\frac{3\,{b}^{3}}{8\,d \left ( a-b \right ) ^{3}{a}^{2}}\arctan \left ({b\cot \left ( dx+c \right ){\frac{1}{\sqrt{ab}}}} \right ){\frac{1}{\sqrt{ab}}}} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

int(1/(a+b*cot(d*x+c)^2)^3,x)

[Out]

-1/2/d/(a-b)^3*Pi+1/d/(a-b)^3*arccot(cot(d*x+c))+7/8/d*b^2/(a-b)^3/(a+b*cot(d*x+c)^2)^2*cot(d*x+c)^3-5/4/d*b^3
/(a-b)^3/(a+b*cot(d*x+c)^2)^2/a*cot(d*x+c)^3+3/8/d*b^4/(a-b)^3/(a+b*cot(d*x+c)^2)^2/a^2*cot(d*x+c)^3+9/8/d*b/(
a-b)^3/(a+b*cot(d*x+c)^2)^2*a*cot(d*x+c)-7/4/d*b^2/(a-b)^3/(a+b*cot(d*x+c)^2)^2*cot(d*x+c)+5/8/d*b^3/(a-b)^3/(
a+b*cot(d*x+c)^2)^2/a*cot(d*x+c)+15/8/d*b/(a-b)^3/(a*b)^(1/2)*arctan(cot(d*x+c)*b/(a*b)^(1/2))-5/4/d*b^2/(a-b)
^3/a/(a*b)^(1/2)*arctan(cot(d*x+c)*b/(a*b)^(1/2))+3/8/d*b^3/(a-b)^3/a^2/(a*b)^(1/2)*arctan(cot(d*x+c)*b/(a*b)^
(1/2))

________________________________________________________________________________________

Maxima [F(-2)]  time = 0., size = 0, normalized size = 0. \begin{align*} \text{Exception raised: ValueError} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(1/(a+b*cot(d*x+c)^2)^3,x, algorithm="maxima")

[Out]

Exception raised: ValueError

________________________________________________________________________________________

Fricas [B]  time = 2.25495, size = 2361, normalized size = 15.74 \begin{align*} \text{result too large to display} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(1/(a+b*cot(d*x+c)^2)^3,x, algorithm="fricas")

[Out]

[1/32*(32*(a^4 - 2*a^3*b + a^2*b^2)*d*x*cos(2*d*x + 2*c)^2 - 64*(a^4 - a^2*b^2)*d*x*cos(2*d*x + 2*c) + 32*(a^4
 + 2*a^3*b + a^2*b^2)*d*x - (15*a^4 + 20*a^3*b - 2*a^2*b^2 - 4*a*b^3 + 3*b^4 + (15*a^4 - 40*a^3*b + 38*a^2*b^2
 - 16*a*b^3 + 3*b^4)*cos(2*d*x + 2*c)^2 - 2*(15*a^4 - 10*a^3*b - 12*a^2*b^2 + 10*a*b^3 - 3*b^4)*cos(2*d*x + 2*
c))*sqrt(-b/a)*log(((a^2 + 6*a*b + b^2)*cos(2*d*x + 2*c)^2 + 4*(a^2 - a*b - (a^2 + a*b)*cos(2*d*x + 2*c))*sqrt
(-b/a)*sin(2*d*x + 2*c) + a^2 - 6*a*b + b^2 - 2*(a^2 - b^2)*cos(2*d*x + 2*c))/((a^2 - 2*a*b + b^2)*cos(2*d*x +
 2*c)^2 + a^2 + 2*a*b + b^2 - 2*(a^2 - b^2)*cos(2*d*x + 2*c))) + 4*(9*a^3*b - 7*a^2*b^2 - 5*a*b^3 + 3*b^4 - 3*
(3*a^3*b - 7*a^2*b^2 + 5*a*b^3 - b^4)*cos(2*d*x + 2*c))*sin(2*d*x + 2*c))/((a^7 - 5*a^6*b + 10*a^5*b^2 - 10*a^
4*b^3 + 5*a^3*b^4 - a^2*b^5)*d*cos(2*d*x + 2*c)^2 - 2*(a^7 - 3*a^6*b + 2*a^5*b^2 + 2*a^4*b^3 - 3*a^3*b^4 + a^2
*b^5)*d*cos(2*d*x + 2*c) + (a^7 - a^6*b - 2*a^5*b^2 + 2*a^4*b^3 + a^3*b^4 - a^2*b^5)*d), 1/16*(16*(a^4 - 2*a^3
*b + a^2*b^2)*d*x*cos(2*d*x + 2*c)^2 - 32*(a^4 - a^2*b^2)*d*x*cos(2*d*x + 2*c) + 16*(a^4 + 2*a^3*b + a^2*b^2)*
d*x + (15*a^4 + 20*a^3*b - 2*a^2*b^2 - 4*a*b^3 + 3*b^4 + (15*a^4 - 40*a^3*b + 38*a^2*b^2 - 16*a*b^3 + 3*b^4)*c
os(2*d*x + 2*c)^2 - 2*(15*a^4 - 10*a^3*b - 12*a^2*b^2 + 10*a*b^3 - 3*b^4)*cos(2*d*x + 2*c))*sqrt(b/a)*arctan(1
/2*((a + b)*cos(2*d*x + 2*c) - a + b)*sqrt(b/a)/(b*sin(2*d*x + 2*c))) + 2*(9*a^3*b - 7*a^2*b^2 - 5*a*b^3 + 3*b
^4 - 3*(3*a^3*b - 7*a^2*b^2 + 5*a*b^3 - b^4)*cos(2*d*x + 2*c))*sin(2*d*x + 2*c))/((a^7 - 5*a^6*b + 10*a^5*b^2
- 10*a^4*b^3 + 5*a^3*b^4 - a^2*b^5)*d*cos(2*d*x + 2*c)^2 - 2*(a^7 - 3*a^6*b + 2*a^5*b^2 + 2*a^4*b^3 - 3*a^3*b^
4 + a^2*b^5)*d*cos(2*d*x + 2*c) + (a^7 - a^6*b - 2*a^5*b^2 + 2*a^4*b^3 + a^3*b^4 - a^2*b^5)*d)]

________________________________________________________________________________________

Sympy [F(-1)]  time = 0., size = 0, normalized size = 0. \begin{align*} \text{Timed out} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(1/(a+b*cot(d*x+c)**2)**3,x)

[Out]

Timed out

________________________________________________________________________________________

Giac [A]  time = 1.24144, size = 278, normalized size = 1.85 \begin{align*} -\frac{\frac{{\left (15 \, a^{2} b - 10 \, a b^{2} + 3 \, b^{3}\right )}{\left (\pi \left \lfloor \frac{d x + c}{\pi } + \frac{1}{2} \right \rfloor \mathrm{sgn}\left (a\right ) + \arctan \left (\frac{a \tan \left (d x + c\right )}{\sqrt{a b}}\right )\right )}}{{\left (a^{5} - 3 \, a^{4} b + 3 \, a^{3} b^{2} - a^{2} b^{3}\right )} \sqrt{a b}} - \frac{8 \,{\left (d x + c\right )}}{a^{3} - 3 \, a^{2} b + 3 \, a b^{2} - b^{3}} - \frac{9 \, a^{2} b \tan \left (d x + c\right )^{3} - 5 \, a b^{2} \tan \left (d x + c\right )^{3} + 7 \, a b^{2} \tan \left (d x + c\right ) - 3 \, b^{3} \tan \left (d x + c\right )}{{\left (a^{4} - 2 \, a^{3} b + a^{2} b^{2}\right )}{\left (a \tan \left (d x + c\right )^{2} + b\right )}^{2}}}{8 \, d} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(1/(a+b*cot(d*x+c)^2)^3,x, algorithm="giac")

[Out]

-1/8*((15*a^2*b - 10*a*b^2 + 3*b^3)*(pi*floor((d*x + c)/pi + 1/2)*sgn(a) + arctan(a*tan(d*x + c)/sqrt(a*b)))/(
(a^5 - 3*a^4*b + 3*a^3*b^2 - a^2*b^3)*sqrt(a*b)) - 8*(d*x + c)/(a^3 - 3*a^2*b + 3*a*b^2 - b^3) - (9*a^2*b*tan(
d*x + c)^3 - 5*a*b^2*tan(d*x + c)^3 + 7*a*b^2*tan(d*x + c) - 3*b^3*tan(d*x + c))/((a^4 - 2*a^3*b + a^2*b^2)*(a
*tan(d*x + c)^2 + b)^2))/d